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Wildfires contribute significantly to global soot emissions, yet their aerosol formation mechanisms and 
resulting particle properties are poorly understood and parameterized in climate models. The conventional 
view holds that soot is formed via the cluster-dilute aggregation mechanism in wildfires and emitted as 
aggregates with fractal dimension D/~ 1.8 mobility diameter D m < 1 ^im, and aerodynamic diameter D a < 
300 nm. Here we report the ubiquitous presence of soot superaggregates (SAs) in the outflow from a major 
wildfire in India. SAs are porous, low-density aggregates of cluster-dilute aggregates with characteristic 
D/~ 2.6, D m > 1 fim, and D a < 300 nm that form via the cluster-dense aggregation mechanism. We present 
additional observations of soot SAs in wildfire smoke-laden air masses over Northern California, New 
Mexico, and Mexico City. We estimate that SAs contribute, per unit optical depth, up to 35% less 
atmospheric warming than freshly- emitted (£>/• — 1.8) aggregates, and —90% more warming than the 
volume-equivalent spherical soot particles simulated in climate models. 



On a global scale, wildfires emit approximately 34% of total atmospheric soot mass, while in certain 
regions such as southeast Asia and Russia, these fires contribute as much as 63% of regional soot mass 
emissions 1 . In the context of climate change, soot emitted from wildfires and biomass burning episodes 
contribute to one of the largest uncertainties in current estimates of radiative forcing 2 . This large uncertainty is 
due to poor understanding of the microphysical properties of wildfire-emitted soot and their parameterizations in 
models and satellite retrieval algorithms 2,3 . In recent years, researchers have made thorough efforts to characterize 
these properties for soot emitted from small-scale, controlled combustion systems— such as prescribed and slash 
burns — as a function of various process parameters such as fuel type, flame temperature, environmental condi- 
tions driving the combustion process, and interrelationships of these parameters 4 " 7 . These studies have been 
conducted under the assumption that they may closely mimic natural wildfires in their soot formation mech- 
anism and emitted particle properties. 

Soot formation in small-scale combustion systems takes place in the cluster- dilute aggregation regime, wherein 
average aggregate-aggregate separation in the aerosol system is much larger than aggregate size 8 . Aggregation in 
this regime proceeds via a three-dimensional, diffusion -limited growth mechanism, governed by the mean-field 
Smoluchowski equation 8 ' 9 . The resultant morphology of aggregates emitted into the atmosphere from this regime 
is described by a fractal dimension Df ~ 1.8, maximum lengths up to 1 urn, mobility diameter D m < 1 um, and 
aerodynamic diameter D a — used for estimating the probability of deposition within lungs 10 — less than 300 nm 5 . 
These aggregates consist of tens to a few hundred monomers, with each monomer typically 30-50 nm in 
diameter 5 . Through atmospheric processing, it is possible for these aggregates to collapse into sphere-like 
morphologies with Df ~ 2.6 11,12 . Climate modelers assume the morphology of bare soot aggregates to be spherical 
in model parameterizations 213 . 

The unpredictable occurrence and nature of large-scale wildfires severely limit the opportunity for in-plume 
sampling of smoke particles for analysis of particle properties and understanding of soot formation mechan- 
isms 3 ' 4 . Single-particle electron microscopy analysis provides crucial information toward accurate and size- 
unconstrained characterization of aerosol microphysical properties 514 . Additionally, information on the aerosol 
formation mechanism can be deduced from single-particle electron micrographs 8 ' 91516 . There have been only a 
limited number of studies conducted to investigate aerosols emitted from natural wildfires using electron micro- 
scopy 17 " 23 . The majority of these studies sampled particles during a fire's smoldering combustion phase (low- 
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Lat 39.20 33.14 27.06 20.96 14.86 8.75 2.63 -3.48 -9.55 

Lon 82.66 80.89 79.30 77.82 76.43 75.09 73.78 72.48 71.17 



Figure 1 | The 2012 Nagarhole forest wildfire smoke plume transport. (A) Average of NOAA HYSPLIT ensemble trajectories ending at the Maldives 
Climate Observatory-Hanimaadhoo, Maldives (MCOH) on 01, 02, and 03 March 2012. (Image created using Adobe Photoshop™); (B) Visible imagery of 
the Indian Peninsula from the MODIS sensor aboard the Terra satellite for 28 February 2012. (Image obtained from NASA Near Real-time (NRT) data 
archive); (C) 532 nm backscatter return signal from the CALIOP Lidar aboard the CALIPSO satellite showing vertical distribution of aerosols (Image 
obtained from NASA CALIPSO data archive). The color scale on the right indicates the strength of the LIDAR return signal: boundary layer clouds usually 
show up as grey or white; cirrus clouds range from yellow to grey; and aerosols show up as green, yellow, and red (indicating low, medium, and high 
loadings, respectively). 



temperature). Exceptions have been the aircraft sampling studies 
conducted during the 1990's in the over-fire regions of flaming forest 
fires in Brazil and southern Africa 1718 . Observations of unusually 
large soot aggregates were made from these fires, but the investiga- 
tors failed to distinguish the microphysical properties of these part- 
icles from conventional sub-micron soot aggregates. As a result, 
these unique observations have gone unnoticed, and there has been 
no follow-up investigation conducted on the formation mechanism 
and frequency of occurrence of these unusually large soot aggregates, 
their microphysical properties, and their potential impact on radi- 
ative forcing and health. 

Here, we investigate particles contained in the flaming-phase 
plumes of the Nagarhole National Forest fire (NNFF) 24 in 
Karnataka (India) and find the ubiquitous occurrence of superag- 
gregates (SAs), a hitherto unrecognized form of soot distinct from 
conventional sub-micron aggregates. We report additional obser- 
vation of these SAs in wildfire smoke-laden air masses over 
Sacramento (Northern California, USA), Los Alamos (New 
Mexico, USA), and the Mexico City metropolitan area (Mexico), 
respectively. Based on the unique morphological properties of SAs, 
we discuss their possible formation mechanism and their potential 
impact on human health. We also compute numerically- exact optical 
properties of these particles and compare them with those of sub- 
micron size soot particles. We make use of the optical properties to 
calculate direct radiative forcing efficiencies of SAs at the top of the 
atmosphere and discuss their net warming or cooling of the atmo- 
sphere. Finally, we address the need for future research to better 
understand and characterize the detection and atmospheric proces- 
sing of soot SAs for quantitatively estimating their impact on climate 
and health. 



Results 

We collected aerosol samples for scanning electron microscopy (SEM) 
analysis downwind of the NNFF over the Indian Ocean at the Maldives 
Climate Observatory on Hanimaadhoo Island (MCOH) (6.78° N, 
73.18° E). The NNFF, which lasted for a week beginning February 
27, 2012, burned approximately 35 km 2 of dry deciduous forest con- 
taining dry bamboo and teak trees. Dense smoke from intense flaming 
combustion was reported, with the event turning into a firestorm 
within a day 24 . The MCOH aerosol number concentration increased 
from about 800 to 3000 cm" 3 during this period. The months of 
November through May constitute the dry season in South Asia, when 
low-level flow brings a polluted air mass from Asia to the Indian 
Ocean 25 . Ensemble back- trajectory analyses (Fig. la; Fig. SI in supple- 
mentary information) coupled with satellite imagery and the Cloud 
Aerosol Lidar and Infrared Pathfinder Satellite (CALIPSO) measure- 
ments (Fig. 1 b and c) show a low-level polluted air mass— between 1 
to 3 km above sea level— transported from the forest fire site flowing 
southwest over the Indian Ocean. Gas chromatography interfaced with 
mass spectrometry analyses of aerosol samples, revealed the presence 
of levoglucosan — a molecular marker for biomass burning emis- 
sions 26 — in trace amounts (about 0.09 ng/m 3 ). 

To investigate how commonly these SAs occur in different geo- 
graphical locations and atmospheric conditions, we sampled aerosols 
contained in wildfire smoke-laden air masses over Sacramento dur- 
ing the CARES (Carbonaceous Aerosol and Radiative Effects 
Study) 27 in June 2010 and over Mexico City as part of the 
MILAGRO (Megacity Initiative: Local And Global Research 
Observations) 28 ' 29 study during March- June 2006. Finally, in 2011 
at Los Alamos, we sampled the downwind plumes of the Las 
Conchas wildfire 19 , the second largest wildfire in the state's history. 
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Figure 2 | Soot superaggregate morphology and particle formation 
mechanism in wildfires. (A, B) SEM images of typical soot SAs observed 
from sampling of smoke plumes from the Nagarhole forest fire. These 
aggregates were greater than 5 microns in maximum length, made of 
thousands of 45-50 nm diameter carbon monomers, and had a 
characteristic cluster-dense fractal dimension Df~ 2.6. (C) A typical sub- 
micrometer size aggregate, made of a few tens of monomers, with a 
characteristic cluster- dilute 1.8. (D) Magnified micrograph of a soot 
SA from this study, highlighting the low degree of organic carbon coating 
between monomers and showing no visible evidence of atmospheric 
processing during long-range transport. 

Morphological properties of soot superaggregates. We measured 
the structural and fractal properties of individual carbonaceous 
particles collected from smoke plumes over the four locations 
using established image analysis routines 5 . The fractal properties of 
individual particles were quantified using the perimeter and 
ensemble methods 30 " 32 . Electron micrographs of typical SAs 
observed at the four sampling sites are shown in Fig. 2 and 3. The 
SAs consisted of aggregates of sub -micrometer size, cluster- dilute 
soot aggregates with characteristic Df = 1.9 ± 0.2 (Fig. 2c). The 
mean maximum length scales of SAs were between 10 and 20 urn, 
and SAs had distinct D f = 2.6 ± 0.1. 

Our analysis showed that a typical SA consisted of around 3000 
monomers, after accounting for apparent monomer overlap, which is 
parameterized by a power-law factor 1.09 33 . Observation of Df = 1.9 
sub -micron aggregates within individual SAs confirms that the SAs 
were formed via percolation of these aggregates in the fires. High 
magnification images (Fig. 2d) of NNFF revealed minimal coating of 
condensed organic matter on the monomers. This suggested that 
these particles were formed under near unity net equivalence ratio 34 
and resisted atmospheric processing during long-range transport 11 . 
Alternatively, it could be that there was not enough condensable 
organic matter available in the flaming fire plumes of NNFF (see 
table SI in supplementary information) to coat the SAs. 

We analyzed 69 individual carbonaceous particles collected on 
SEM filters from the NNFF. Approximately 99% and 75% of the 
particle mass and number, respectively, were soot SAs 9,31 ' 32 ' 35 , with 
the remaining being aggregates (Fig. 4). We probed the elemental 
composition of the particles using energy dispersive X-ray spectro- 
scopy (EDX), finding carbon and oxygen to be their primary con- 
stituents. No tar balls 22 ' 36 or particles with inorganic inclusions were 
observed. We analyzed 580 particles from the CARES campaign and 
found approximately 16% of the total aerosol particle number was 
SAs. The remaining 57% and 27% of particles were externally sub- 
micron soot (occurring as bare and mixed with organic carbon) and 
dust particles, respectively. The SAs observed in samples from 
MILAGRO and the Las Conchas fires were less than 1% in number. 
Aerosol types in these two locations were mostly organic carbon and 
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Figure 3 | Soot superaggregates observed in biomass burning plumes in 
different locations in USA and Mexico. (A, B) Typical soot SAs observed 
in the Sacramento urban area. The sampled air mass containing biomass 
burning smoke particles is believed to be transported from the neighboring 
Sierra Nevada foothills and San Joaquin County 27 ' 57 (C) The Las Conchas 
fire was the second largest wildfire in the state history of New Mexico 
burning an area of 634 km 2 . The sampled smoke plumes from this fire for 
this study were mostly from smoldering phase mixed with intermittent 
naming phase. As a result, we observed majority (—80%) of particles to be 
spherical tar balls, observable in the background of the soot SA particle 19 . 
(D) A soot SA sampled in Mexico City. The air mass containing these 
particles is believed to originate from burning of pine savannas in the 
mountains near Mexico City 28 . 

tar balls 19 ' 20,28 ' 29 ' 37 , suggesting dominant emission from the smolder- 
ing phase of wildfires. 

For the SAs and cluster- dilute aggregates observed at MCOH and 
during CARES, we calculated their mobility diameters D m 14,38 and 
mass distributions (Fig. 4 and 5) based on their single-particle pro- 
jected area equivalent diameters 14 . The SAs had a range of D m 
between 1 and 20 |im with a mean D m ~ 3 |im. The monomer 
number size distribution could be described by a mono-modal log- 
normal size distribution with a mean monomer diameter of 50 nm 
and a standard deviation of 5 nm. We observed a majority (95%) of 
SAs in the third stage, D a < 0.3 |im, of the impactor used for aerosol 
collection. Although characterized by very large geometric diameters 
(D m ), the low D a of SAs suggest that they are highly porous, have low 
effective densities, and could get deposited in the innermost lung 
airways and alveoli via the process of diffusion deposition, similar 
to soot aggregates 10,39 . However, the extent of lung penetration of SAs 
depends on particle D a 40 . Super-micron size porous carbon aggre- 
gates, similar to soot SAs in this study, have been synthesized in the 
laboratory 41 , and were observed to have effective particle densities as 
low as 2.5 mg/cm 3 . 

Discussion 

Superaggregate formation mechanism. SAs are formed when 
cluster- dilute aggregates enter into a cluster- dense aggregation 
regime in flames 35 . This regime, defined by a small ratio of the 
mean aggregate nearest-neighbor separation to aggregate size and 
by enhanced kinetics, results in the aggregates sticking together and 
percolating to form a volume spanning SA with a universal Df ~ 
2.6 9 ' 35 . SAs usually consist of more than a few thousand monomers 
and are greater than 1 |im in length 31 ' 32 . Over small length scales, SAs 
have a Df ~ 1.8 because of their formation via percolation of cluster- 
dilute aggregates 9 . 
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Figure 4 | Number distribution of soot aggregates and superaggregates from the Nagarhole forest fire. (A, B) Number size distribution of the particles is 
calculated in terms of the projected area equivalent diameter, which is equivalent to particle mobility diameter. (C, D) Single particle mass 
distribution, calculated from two-dimensional particle images using a carbon monomer density value =1.8 g/cm 3 . It is to be noted that sea-salt particles 
from the Indian Ocean accompanied the carbonaceous aerosols and were excluded from this analysis. 



With the exception of large-scale sooty and turbulent fires 32 , 
upward- rising flame systems cannot emit SAs into the atmosphere 
due to their narrowing flame fronts and associated buoyancy- driven 
aggregate fragmentation effects 31 . For SAs to be emitted as final 
products, an aggregating system must maintain the right physical 
conditions in terms of threshold particle volume fractions and res- 
idence time. A large-scale turbulent flame system provides these 
conditions by trapping sub-micron aggregates in its vortices 9 ' 31,32,42 
(Fig. 6). Once trapped, these aggregates begin to jam together under 
high particle volume fraction conditions to form SAs. For the soot 
SAs observed in this study, we hypothesize that they were formed in 
the vortices of the wildfire flame bodies and subsequently emitted to 
the atmosphere. 

Impact on direct radiative forcing. Quantitative knowledge of the 
soot morphologies enabled numerically- exact superposition T- 
matrix modeling of single -particle optical properties— asymmetry 
parameter g and single scattering albedo co — at 550 nm as 
functions of aggregate Df and monomer number (Fig. 7a and b). 
We also computed these parameters using Mie theory for volume- 
equivalent spherical particles 14 as used by current climate models. 
Aggregate simulations and their T-matrix modeling methods are 
described in detail elsewhere 43 . Note that the T-matrix calculations 
required significant computational power and time, limiting the total 
number of monomers that could be investigated to 1000, but it is 
likely that these results are broadly representative of SA optical 
properties. Both g and co are fundamental parameters for aerosol 
radiative transfer calculations 44 , co is the ratio of particle scattering 
to extinction cross section, and ranges between 0 for a purely 
absorbing particle to +1 for a purely scattering particle. The 
asymmetry parameter (g) is defined as the intensity-weighted 
average cosine of the scattering and ranges between — 1 for entirely 



backscattered light to +1 for entirely forward scattered light, co 
increases with increasing monomer number as well as with Df in 
both the aggregate and SA regimes. For Mie-spheres, their co 
increases with increasing size. SAs have higher co values than 
freshly emitted aggregates (Df = 1.8) and slightly lower and higher 
values than aged aggregates (Df ~ 3) for monomer numbers 400-700 
and >700, respectively. Alternatively, SAs have higher g values than 
both aggregates and Mie-spheres. The g values for aggregates 
increase monotonically and reach a maximum value at Df ~ 2.75, 
after which they start to slowly decline to settle at a lower value for 
Df = 3.0. Mie-spheres show even lower g values than Df = 1.8 
aggregates. This behavior can be attributed to an increase in 
backscattering by particles as they approach near-spherical 
morphology 44 . 

Using these calculated single-particle optical properties, we esti- 
mated the top of the atmosphere direct forcing efficiency (DFE; the 
radiative forcing per unit optical depth and bandwidth 45 ) at 550 nm 
for aggregates, SAs, and Mie-spheres as functions of Df and size 
(Fig. 7c). We chose the complex refractive index of soot to be 
1.95-0.79i and the up-scatter fraction as a function of g per past 
recommendations 46,47 . The difference in the DFE value between 
freshly emitted (i.e. Df = 1.8) 11 aggregates and SAs is up to 
+ 7 Wm" 2 . This positive difference in forcing is the excessive short- 
wave heating of an optically thin layer of the troposphere by freshly 
emitted cluster- dilute aggregates compared to SAs. The forcing dif- 
ference between SAs and atmospherically processed or aged aggre- 
gates (Df = 2.5-2.7) 11,12 is negligible. However, in comparison with 
completely collapsed 48 (Df ~ 3.0) aggregates, SAs contribute to an 
enhanced warming effect by —25% (+4 Wm" 2 ). Note that a Df~ 3.0 
aggregate shouldn't be confused with a volume- equivalent Mie 
sphere, as simulated in climate models. The morphology of a 
Df = 3.0 aggregate differs from a volume-equivalent sphere. A sphere 
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Figure 5 | Number distribution of soot aggregates and superaggregates observed during CARES. The projected area equivalent diameter number 
distribution (A, B) and mass distribution (C, D) of soot particles emitted from forest fires in the Sierra Nevada Foothills and the San Joaquin County 
regions of California, USA. The plumes were sampled in the month of June in an urban location (Sacramento), and also included inorganic aerosols in 
pure and externally- mixed states. These particles were not included in this statistical analysis. 



has a smooth annular region compared to a Df = 3.0 aggregate, which 
has an irregular boundary composed of monomers with voids in 
between them. This difference in morphology significantly affects 
their optical properties 49 . Our calculations show that SAs contribute 
90% more warming compared to volume-equivalent Mie-spheres. 

Future research direction. Taking into account both the abundance 
of soot SAs in the atmosphere and their microphysical properties, 




Figure 6 | Schematic (drawn using Adobe Photoshop™) hypothesizing 
the soot superaggregate formation mechanism in wildfires. Particle 
residence times are enhanced for cluster-dilute aggregates that get trapped 
into one of the many vortices inside a large-scale, turbulent flame. This 
facilitates increase in particle volume fractions and percolation of 
aggregates into superaggregates inside the flame. 



our findings bring to the forefront the significant impacts of these 
previously unrecognized pollutants on climate, human health, air 
pollution monitoring, and mitigation strategies. The higher heating 
effect of these particles, compared to volume- equivalent spheres, 
could change current estimates of climate forcing by models. 
Observation of these particles also demands new investigations on 
multiple research fronts. Instrument development for real-time 
detection and measurement of soot SAs using different size 
conventions, such as mobility, aerodynamic, and volume- equiva- 
lent diameters, is needed. The widely used particle-size monitoring 
instrument, the scanning mobility particle spectrometer 414 , is 
incapable of detecting these particles. The effect of humidity and 
atmospheric processing on soot SAs is currently unknown. Both 
parameters have been shown to affect radiative properties of 
particles 11 . Based on future research findings, effective measures — 
such as appropriate mechanical filtration systems— could be 
employed during large-scale wildfires to help control and mitigate 
public health impacts of soot SAs. 

Methods 

Aerosol sampling. In the Maldives, air was sampled from inlets on a 15 m tall tower 
equipped with various sensors. Aerosol particles were collected using a three-stage 
aerodynamic impactor (MPS-3; California Measurements, Inc.) for single-particle 
SEM analysis. The nominal D a ranges for the three stages were >2, 2-0.3, and 
<0.3 um. Sampling duration ranged from 45 minutes to two hours twice a day, 
during morning and afternoon. Particles were deposited directly onto 10 um thick 
nuclepore clear polycarbonate filters (100-nm diameter pores) applied to metal stubs 
and located at each impactor stage 31 . Aerosol sampling during the CARES, 
MILAGRO, and the Las Conchas wildfire were also conducted on similar nuclepore 
filters using a filtration technique. After exposure, the filters were stored in a dark 
refrigeration unit to avoid particle transformation due to aging. Sampling duration 
for the CARES study was from 2331 (06/21/2010) to 0601 hours (06/22/2010) local 
time, while the sampling duration for the MILAGRO study was conducted from 1359 
to 1546 hours on 03/28/2006. The Las Conchas fire plumes were sampled at Los 
Alamos between 1300 and 1800 hours on 07/22/2011. 
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Figure 7 | Radiative properties of soot superaggregates and comparison 
with aggregates and spherical particles. Optical properties-single 
scattering albedo (A) and asymmetry parameter (B)-and direct forcing 
efficiency (C) dependence of aggregates and superaggregates as functions 
of size (number of monomers) and fractal dimension Df Shaded regions 
represent the upper and lower bound values for Df in the range between 1.8 
and 3.0. The calculations used mean refractive index = 1.95-0.79i, 
monomer diameter = 50 nm, and wavelength = 550 nm. Monomer 
number = 400 corresponds to mobility diameter = 1 (am. 

Microscopy analysis. The filter samples were prepared for SEM analysis by coating 
them with a 1-nm thick layer of platinum to prevent aerosol charging during SEM 
analysis. A field- emission SEM (Hitachi S-4700) was used to analyze the coated filters 
for individual particle morphology and EDX analysis. A relatively moderate 
accelerating voltage of 20 kV was used for imaging most particles. Past studies have 
shown that imaging at this operating voltage has negligible impact on aggregate shape 
due to charging 5,33 . Images of individual particles were selected, based on the random 
distribution of particles on the SEM filter, for shape and size quantification. 

The perimeter fractal analysis method. Calculating single-particle fractal dimension 
with the Perimeter Method is done by drawing grids of differing box sizes upon a 2-d 



image of a fractal aggregate. Given a grid with a certain box size, the number of 
grid boxes through which the perimeter passes is counted. A grid with a different box 
size is then drawn, and the number of grid boxes through which the perimeter passes 
is counted once again. This process is repeated, and the logarithm of the box size 
plotted against the logarithm of the box count, providing the parameter fractal 
dimension D p as the slope. Past experimental and simulation studies 3130,51 have 
confirmed the D p of a superaggregate to be well defined and related to its mass fractal 
dimension Dousing the following empirical relationship: 



D p = l+ 3-D, 



,3/2 



D f >2, 



D 0 = Df, D f <2. 



(1) 



(2) 



The ensemble fractal analysis method. This method is the most accurate method for 
calculating Df of an aggregate. Calculating aggregate with this method involves 
using the relationship expressed by Equation 3 below and determining the values for 
the number of monomers AT with the knowledge of the projected area of the aggregate 
(A a ), the mean projected area of the monomers (A p ), prefactor k 0 (typically 
approximated by a value ~ 1), and a = 1.09 



N = k Q 



(3) 



Aggregate generation and T-matrix calculation. The aggregate simulation method 
involved generating a 3-d fractal aggregate by the particle- cluster aggregation 
technique 52 with a sequential algorithm that intrinsically satisfies the fundamental 
fractal equation 53 . For pre-specified values of and prefactor (typically between 1.1 
and 1.3) 54 , the aggregate generation process is initiated by randomly attaching two 
monomers to each other, followed by the controlled addition of further monomers to 
the cluster at specific positions fulfilling the following conditions: 1) the monomers do 
not intersect, i.e., they make point contact, and 2) the radius of gyration of the new 
aggregate (calculated based on the known positions of the monomers) satisfies the 
fundamental fractal equation for the fractal dimension and pre-factor selected. In this 
study, the pre-factor was set at 1.18 for all fractals as recommended by Oh and 
Sorensen 54 , and the monomer diameter was set to 50 nm for all calculations, based on 
electron microscopy observation of superaggregates in this study. We calculate the 
optical properties of fractal-like soot aggregates using the efficient superposition T- 
matrix code developed for multi- sphere groups with random orientation. The code is 
documented 55 and is available on the Web, free of charge 56 . The critical advantages of 
this method are that it is numerically exact and is much more efficient than any other 
numerical technique based on an explicit solution of the Maxwell equations. 

Aerosol direct forcing efficiency calculation. The shortwave aerosol DFE at the top 

of the atmosphere caused by a uniform, optically thin aerosol layer in the lower 
troposphere was calculated using 



AF 



-SD(l-A dd )T 2 atm (l-R sfc ) 2 



2R sfc 



I — CO 



-001 



(4) 



where AF is the change in net solar flux at the top of the atmosphere due to the 
presence of the aerosols, t is the aerosol optical depth, S is the solar constant, set to 
1370 Wm" 2 , D is the fractional day length, set to 0.5, A dd is the fractional cloud cover, 
set to 0.6, T atm is the solar atmospheric transmittance, set to 0.76, R s f c is the surface 
albedo, set to 0.15 (appropriate for an urban area), and co is the aerosol single 
scattering albedo. The parameter /? is the up-scatter fraction, which is a function of 
asymmetry parameter g as follows: 



= 0.50 - 0.45g + 0.25g z - 0.29^ 



In all of the calculations, the relative humidity is assumed to be 0%. 



(5) 
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